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Abstract 

A new approach to obtain the closure formulas for the von Karman-Howarth 
and Corrsin equations is presented, which is based on the Lagrangian repre¬ 
sentation of the fluid motion, and on the Liouville theorem associated to the 
kinematics of a pair of fluid particles. This kinematics is characterized by the 
hnite-scale separation vector which is assumed to be statistically independent 
from the velocity held. Such assumption is justihed by the hypothesis of fully 
developed turbulence and by the property that this vector varies much more 
rapidly than the velocity held. This formulation leads to the closure formu¬ 
las of von Karman-Howarth and Corrsin equations in terms of longitudinal 
velocity and temperature correlations following a demonstration completely 
diherent with respect to the previous works. Some of the properties and 
the limitations of the closed equations are discussed. In particular, we show 
that the times of evolution of the developed kinetic energy and temperature 
spectra are hnite quantities which depend on the initial conditions. 
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1. Introduction 

Recently, a work dealing with the bifurcations analysis of the turbulent 
energy cascade llj] presents, among the other things, the relationship be¬ 
tween Navier-Stokes equations bifurcations and turbulent energy cascade, 
showing that these bifurcations produce a negative value of the skewness of 
the longitudinal velocity difference, and a separation rate between contiguous 
trajectories which exponentially diverges with the time. 

The present work proposes a specific analysis of the isotropic homoge¬ 
neous turbulence of incompressible fluids based on Ref. jll|, and on the 
property that the contiguous fluid particles trajectories continuously diverge 
due to the Navier-Stokes bifurcations. The proposed formulation adopts the 
Lagrangian representation of the fluid motion and the Liouville theorem. 

According to the present study, the bifurcations determine the energy 
cascade, where the separation vector ^ between two fluid particles trajecto¬ 
ries varies much more rapidly than the velocity and temperature ll|. Such 
property, in conjunction with the hypothesis of fully developed turbulence, 
justifies the assumption that ^ and the velocity held are statistically inde¬ 
pendent. This is the crucial hypothesis of the present analysis that allows, 
through the Liouville theorem, to analytically express the closure formulas 
of the von Karman-Howarth and Corrsin equations in terms of longitudinal 
velocity and temperature correlations. The proposed analysis also quantihes 
the mean effect of the trajectories separation through the average hnite- 
scale Lyapunov exponent, a quantity depending on r which gives the rate of 
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separation of the trajectories with hnite distance. This exponent is here cal¬ 
culated in function of the maximum hnite-scale Lyapunov exponent through 
the distribution function of ^ and the Liouville theorem. 

For sake of the reader convenience, we report the main keypoints of the 
article: The hrst part of the work, devoted to the Lagrangian description 
of the fluid motion, provides the representation of velocity and temperature 
helds, and of the kinematics of a pair of fluid particles. The Liouville theorem 
is then introduced, and the statistical independence of ^ from velocity and 
temperature helds is properly justihed. In the second part, a relationship 
between average and maximal Lyapunov exponents of hnite-scale, useful for 
the subsequent analysis, is determined, and the average Lyapunov exponent 
is expressed in terms of longitudinal velocity correlation. Thereafter, the 
closure equations are obtained through the analytical elements introduced in 
the previous two parts. 

The obtained results agree with those presented in Refs. Q Q] , where 
the closure formulas are carried out using a fully diherent formulation and 
exploiting some of the properties of motion of the hnite-scale Lyapunov ba¬ 
sis, and the frame invariance of the triple correlations appearing in the von 
Karman-Howarth and Corrsin equations. This corroborates the results ob¬ 


tained in Refs. 


12 


, ll3|, showing t 


proach and the analysis of Refs. 


;he equivalence between the present ap- 
12|,ll3|. 

Finally, some of the properties of the closed von Karman-Howarth and 
Corrsin equations are studied and their limits of validity are discussed. Specif¬ 
ically, we show that the times of developing of velocity and temperature 
correlations are both hnite quantities which depend on the initial conditions. 
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2. Background: longitudinal velocity and temperature correla¬ 


tions 


In fully developed homogeneous isotropic turbulence, the turbulent fluc¬ 
tuations of velocity and temperature are represented by the following pair 

ibi 

{UrU'j.) 


correlation functions 


20 


/ = 


fo = 




( 1 ) 


(M) 


where the brackets {) denote the average calculated over the ensemble of 
velocity and temperature in the points x and x -|- r, = u(x) • r/r and 
u'^ = u(x-|-r) - r/r are the longitudinal components of fluid velocities, ■d and 
■d' are the fluid temperatures in x and x -|- r, and = {u ■ u)/3, = (d^) 


represent the standard deviations o 
in space due to homogeneity [20|, l9| . Such corre 
von Karman-Howarth and Corrsin equations 


velocity and temperature, both constants 
ations vary according to the 


3, y, y, 


lOj which are the 


evolution equations for / and fo, respectively. These equations, obtained 


throng 


are 
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where v and y = k/{pCp) are kinematic viscosity and thermal diffusivity, 
Cp and k are specihc heat at constant pressure and thermal conductivity. 
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respectively. The quantities Xt = \/—1//"(0) and Xq = 2//^ (0) are the 

Taylor and Corrsin microscales. The boundary conditions associated to Eqs. 
(Ej) and ([3]) are 


/(O) = 1, lim /(r) = 0, 

r—)-oo 

(4) 

fg{0) = 1, lim fg{r) = 0, 
r—>-cxD 


In Eqs. ([2]) and ([3]), K and G arise from the inertia forces and from the 


convective terms of the Navier-Stokes and thermal energy equations 
and can be expressed as 

d 


20 


a, 


K = 


G = 


dvk 

A 

dri. 


{UrU'^ (Ufc -Uk)) , 

{m' (u'fc -Uk)), 


(5) 


Accordingly, K and G do not modify neither the kinetic energy nor the 
thermal energy, and provide the energy cascade mechanism the effect of which 
vanishes for r = 0 and for r —)■ oo, i.e. 


A'(O) = 0, lim K{r) = 0, 

r^oo 

( 6 ) 


G(0) = 0, lim G{r) = 0 

r—>-oo 


In this study, we analyze the case where u does not depend on t?, therefore 
Eq. (EJ is independent from Eq. ([3]) and fg. Conversely, the temperature 
fluctuations will depend on u. 
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3. Lagrangian representation of flnid motion, Lionville theorem 


Following the present analysis, the turbulence is caused by the bifurca¬ 
tions of the Navier-Stokes equations, whereas the temperature plays the role 
of the passive scalar. This study applies also to any passive scalar which 
exhibits diffusivity. These bifurcations frequentely occur in developed tur¬ 
bulence determining a condition of fully developed chaos where velocity and 
temperature exhibit chaotic fluctuations, and the contiguous fluid particles 
trajectories diverge continuously with exponential growth rate ll(|. This 
implies that the fluid deformations are represented by exponential growth 
functions of the time, whereas the fluid state variables, like velocity and 
temperature, are slow growth functions of t. 


3.1. Representation of velocity and temperature fields 

To represent velocity and temperature fluctuations, we start from Ref. 
jn|, where the Navier-Stokes equations are written in the symbolic form 
of operators. Here, the Navier-Stokes equations are expressed following the 
Lagrangian description of motion 32|, and to analyze the temperature fluc¬ 
tuations, the thermal energy equation is also given in the same form 


u = N(u) = Np(u) -I- z/Lnu, 


R = M(i?) = xLM'd 


(7) 


The first equation is the Navier-Stokes equations in the symbolic form of 


operators 


33| . where the pressure held p has been eliminated through the 


continuity equation, u and u denote, respectively, velocity and acceleration 
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fields 32|, Ln(o) = Lm(o) = V^(o), and Np(u) = —Vp/p. Specifically, 
Np(u) is the integral non-linear operator which expresses the pressure gra¬ 
dient through the velocity held in the entire huid domain. That is, the 
pressure gives the non-local effect of the velocity held 33|], and the Navier- 


Stokes equations are reduced to be an integro-diherential equation formally 
expressed by Eq. 0 in the symbolic form of operators. 

The other equation describes the huid temperature variations and is also 
in the symbolic form of operators. The temperature held i? is denoted in 
bold type as it is an element of the vector space of the temperature helds 
{■j?} which is an inhnite dimensional manifold, thus d and d' are the values of 
■j? calculated in x and x-|-r, respectively, and 'd is the material time derivative 
of temperature [32 1. 


Now, the bifurcations are responsible for the huctuations of u and whose 
statistics is formally described by the distribution F=F{u, i?}, a functional 
of u and which satishes the Liouville theorem associated to Eqs. m 
OF 

— + Vu • (NF) + V,, • (MF) = 0 (8) 

where the operators Vu ■ (o) and V,? • (o) denote the divergence in the vector 
spaces of velocity and temperature helds. This theorem is derived from Eq. 
()7[) and from the condition that the integral of F over V identically equals 
the unity 


F dV = 1, Vt > 0. 


(9) 


where V = {u} x {"i?} is the phase space associated to Eq. ([7]), in which 
{u} and {"i?} are the vector spaces of velocity and temperature helds, and 
dV = du d'd represents the elemental volume in the phase space V. 
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3.2. Representation of fluid kinematics. Kinematics of a pair of fluid parti¬ 
cles 

To complete the Lagrangian description of fluid motion, the kinematics 
equations are added to Eqs. ([7]). To our purposes, we consider the trajectories 
of fluid particles passing through x and x + r at f =0, whose evolution 
equations are 

X = u(a;, t), 


i = u{x + i,t) - u(x,t), 


( 10 ) 


with the initial conditions 


a;(0) = X, 


( 11 ) 


1(0) = r 

where £ is the separation vector and ^(0) = r is its initial value. This 
separation vector is a function of t of exponential growth llj] which depends 
on the scale r = |r| and on the initial condition ^(0), whereas u varies 
according to Eq. ([7]). Following Eq. ffTH . the map x ■ ^(0) ^ ^(^) 
provides the current separation distance £ between two fluid particles which 
are located at the referential positions x and x + r at t = 0 32| . It is worth 
to remark that x and x + r are hxed points of the inertial frame, whereas x 
and a; + £ represent the positions of the fluid particles which vary with the 
time and that satisfy the initial conditions flTT]) . 

To analyze the average effect of the trajectories separation, the statis¬ 
tics of ^ is now introduced. This statistics is described by the distribution 





function of ^ which does not depend on x due to turbulence homogeneity 


p = p{t,i) 


( 12 ) 


where 



1, Vt > 0 


(13) 


in which H = {a;} x and dS = dx d^ represents the relative elemen¬ 
tal volume. As the problem fllUl) is studied in an inhnite fluid domain, P 
identically vanishes on the boundaries of H, i.e. 


-P(^)^) = 0 y{x,^)EdS, Vt > 0 (14) 


Next, this PDF satishes the Liouville theorem associated to Eqs. flTU]) and 

ca 

dP d / • \ d 

where the differential operators d/d$^ ■ (o) and d/dx ■ (o) indicate the di¬ 
vergence of (o) with respect to ^ and x, respectively. Taking into account 
the hypothesis of homogeneity and the continuity equation, the Liouville 
equation reads as 


dP 



(16) 


As the initial condition of ^ and x is given by Eq. (na, the initial condition 
of the distribution function P = P(t, is 


F(0,O = ^(l,r) 


where (5(.,.) is the Dirac’s delta. 


(17) 
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Therefore, the average of any quantity depending on x and say "0 = 
will be calculated as integral of P over {H}, according to 




P dH 


(18) 


3.3. Statistical independence of ^ from u and 

At this stage of the analysis, it is worth to remark that directly 
responsible for the relative trajectories mixing 28|, behaves like a function 
of fast growth of t, therefore it varies much more rapidly than u and i? 11|. 
Accordingly, in fully developed turbulence, the time-scales of ^ are expected 
to be completely separated from those associated to u and "i? in the sense that 
^ and (u, "i?) exhibit chaotic behavior and their power spectra are supposed 
to be located in frequency intervals which are completely separated. For this 
reason, £ and (u, -i?) are considered to be statistically uncorrelated. From the 
physical point of view, this means that the effect of the trajectories mixing is 
much more rapid and statistically uncorrelated with respect to the dynamics 
of the fluW system. This property is supported by the arguments presented 
(and references therein), where the author observes the that: a) 


in Ref. 


The helds u(a;,f), (and therefore also u(a:: + ^,t) — u{x,t)) produce chaotic 
trajectories also for relatively simple mathematical structure of the right- 
hand sides u(a;,f) (also for steady helds!). b) The hows given by u{x,t) 
(and therefore also u(a: + ^,t) — u{x,t)) stretch and fold continuously and 
rapidly causing an ehective mixing of the particles trajectories. Hence, it 
is reasonable to postulate that ^ and (u, -d) are statistically uncorrelated, 
and that all the quantities depending on ^ are fast variables whose average 
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is calculated through Eq. flTS]) . This is the fundamental hypothesis of the 
present work following which the distribution function of u and is the 
product of P and F\ 

= (19) 


The time derivative of P is then 

dP dP 


dP 

P-\ - F 

dt dt dt 


( 20 ) 


The hrst term at the R.H.S. of Eq. 020 p is due to the time variations of u 
and 'j 9 following Eqs. ([7]), thus such term gives the time variations of kinetic 
and thermal energies. As far as the second term is concerned, it is related 
only to the fluid kinematics in line with Eq. flTBp . thus it does not contribute 
to kinetic and thermal energies variations. 

We conclude this section by observing the limits of validity of Eq. ([19]). 
As these limits arise from the hypothesis of fully developed turbulence, in 
the cases of intermediate stages of turbulence or in decaying turbulence, the 
statistical independence of ^ and u could be not verihed. 


4. Average finite—scale Lyapunov exponent in terms of the maxi¬ 
mal finite—scale Lyapunov exponent 

The purpose of this section is to present the link between average and 
maximum Lyapunov exponents of hnite-scales that will be useful later on in 
this analysis. 

The average effect of the trajectories divergence is quantihed by the av¬ 
erage hnite-scale Lyapunov exponent, a quantity which expresses the sep¬ 
aration rate between two trajectories with hnite distance. To dehne this 
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exponent, observe that the following quantity 

= ( 21 ) 

gives the local growth rate of |^|. Following its definition, the maximal hnite- 
scale Lyapunov exponent gives the maximum growth rate of the trajectories 
distance. In the present framework, this exponent can be written as 

A = (22) 

The over line o denotes the average calculated according to Eq. (ITSH . where 
A = A(r) thanks to the turbulence isotropy. This exponent is associated to 
a given spatial direction /i, variable with the time, which corresponds to the 
maximum growth rate of trajectories distance. 



Figure 1: Scheme of definition of A', A and A'. 

On the other hand, the trajectories separation happens with a mean 
growth rate. This latter is defined through the average finite-scale Lyapunov 
exponent, a quantity calculated as the mean of A' 














If the ergodic hypothesis is satished, this average exponent is calculated in 
terms of the distribution P as 

^ /s^ H = 5 

With reference to Fig. [H the continuous line represents the time variations 
of l/21n(£ ■ and its local slope gives A' (dash-dotted line), whereas A and 
A' are represented by dashed and dotted lines, respectively. 

The two exponents are related with each other, and to express A' in 
function of A, we write A' in terms of dP/dt 

where the first term at the R.H.S. of Eq. (I25ll identifies A, whereas the second 
one is expressed by means of the Liouville theorem 

(26) 

Integrating by parts the last term at the R.H.S. of Eq. fl26l) . and using the 
Green’s second identity, we obtain the sum of two terms. The hrst one of 
these identically vanishes as it is an integral of P over d {^} where P = 0, 
whereas the other one is different from zero. This leads to 


A' = A- 

Hence, taking into account Eq. 


, we have 


(27) 




(28) 


Due to homogeneity, the separation rate does not depend on x, and thanks 
to the isotropy A' and A are both even functions of r 0,Q. 
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Finally, it is worth remarking a property of the function x • £(0) $'■ 

when 0 < t < 1/A, ^ rexp(At), thus x satishes the following equation 

5 ; = ;. 0<t<l/\ (29) 


5. Average Lyapunov exponent in terms of velocity correlation 

In order to express A' in terms of /, observe that along the direction li, 
^ varies following ^ = ,^(0) exp(Af), i.e. 

^ = Q(f)^(0)exp(At) (30) 

where Q(f) is a proper rotation matrix which gives the orientation of ^ with 
respect to the inertial reference frame. Hence, the relative velocity ^ is 

^ X ^ + A|, (31) 

where (jJ^ is the angular velocity of ^ and li with respect to the inertial 
reference. 

Now, because of the isotropy, the standard deviation of Aur does not 
depend on r/r, and this implies that 

(« - Urf) = ■ 0 = AV^, (32) 

Hence, / and A' are linked with each other according to 

Equation fl33|) is defined only for / < 1 and this agrees with Schwarz inequal¬ 
ity j^, and in particular, for df /dr <0 near the origin. 
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In conclusion, in fully developed homogeneous isotropic turbulence, the 
fluid particles trajectories continuously diverge with an average separation 
rate R which depends on r following 

R = J'r = (34) 


6. Closure equations 


Here, the closure formulas of Eqs. (j2]) and ([3]) are determined by means 
of the elements introduced in the previous sections. 

To obtain these formulas, observe that K and G are responsible only 
for the energy cascade. This latter is the kinetic and thermal energy flow 
between the length scales, without changing the total amount of kinetic and 
thermal energy (see for instance j^, 9| and references therein). In fact, inertia 
and pressure forces and the convective term produce an interaction between 
the various Fourier components of velocity and temperature spectra which 
gives the kinetic and thermal energy transfer between the volume elements 
in the wave-number ^ace, where the global effect of such interaction leaves 


and 6*^ unaltered 


, 3 - 


Now, to determine K and G, we first calculate the pair correlations / and 
fe through the PDF R, taking into account the initial condition of P given 
by Eq. (ITT)) . Accordingly, longitudinal velocity and temperature correlations 
can be expressed through R 

M^/(r) = {uru'^) = j R u^{x)u^{x + ^) dVdH = ^ F Uru'^ dV, 

(35) 


n Jv 


e‘^fg(r) = i'dR) = f [ R ^{x)^{x + I) dVdS = [ F M dV, 

Jv Jv 
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or in terms of F as 


wV(0 = + = / F u^{x)u^{x + ^) dV, 

Jv 


(36) 


0'^h(i) = + 0) ^ [ F'd(x)'S(x + i)dV, 

Jv 

in which = u • ^/,^. The correlations fl5^ and fl5B]) are linked with each 
other because the fluctuation x ■ ^(J~) ^{t~) = r, t = t~ + e, e > 0 

varying much more rapidly than / and fg, leaves unaltered the correlations 
in the transformed points, therefore, taking into account Eq. 0291) . we have 
fe{t,0=fe{'t~,r), and 


m,- 

mi 


dfjr] 

dr 


dfgjr] 

dr 


dr 
t- ^ 


dr 


df{r] 

dr 


f dfejr) 
\ dr 


r 


(37) 


r 

e 


Hence, the time derivatives of the correlation functions are calculated by 
means of dF/dt, using Eq. (]20|) 


m = /s X (f ^^«) 


(38) 


As previously seen, the terms with dF/dt correspond to the time variations of 
kinetic and thermal energies and to the changing of the correlations caused 
by viscosity and thermal conductivity, whereas those with dP/dt, arising 
from the fluctuations of do not modify neither vJ nor 6*^. This identifies 
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K and G as 


K = J J F u^{x)u^{x + dVdS 

dP 


(39) 


G = 



n JV 


dt 


F 'd{x)'d{x + dVdS 


where dP/dt is expressed through the Liouville theorem flThl) 

K = - y ^ ^ u^{x)u^{x + dVdH 


G = - y) y ^ ■ (Pi) F «(x)-d(x + i) dVdS 


(40) 


Next, integrating by parts the R.H.S. of Eq. fHOj) . and using the Green’s 
second identity, we obtain that both K and G are the sum of two terms. The 
hrst ones identically vanish as these are integrals of P over d {^} in which 
P = 0, whereas the other ones are different from zero. This leads to 

d 


K = J^P ^ - — {u^{x)u^{x + ^)) dH 

Taking into account the turbulence isotropy, and Eq. (I37)) . we obtain 


^ + ^)) = ^ + ^)) | = ^ {'d^d') ^ 


(41) 


(42) 


It is worth remarking that depending on x and is considered to be a 
fast variable, thus its contribution is calculated with the statistics of P. This 
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leads to express K and G in terms of df /dr and dfe/dr 


K = u r 


P dH, 


dr J's i 


( 43 ) 


G = 0 


2^/e 


P dS 


dr ^ • I 

Hence, in view of Eqs. fIMl) and fl33ll . we obtain the following expressions 


K = u-^ 


G = e^u 


1 ^ 

dr ’ 


(44) 


1 -/ ^ 
2 dr 


Following Eqs. fl44)) . the phenomenon of kinetic and thermal energy cascade 
is due to the combined action of the spatial variations of /, fe and of the 
trajectories separation. Therefore, K and G result to be non-monotonic 
even functions of r which vanish for r = 0 and tend to zero as r —)■ cx). In 
particular, K gives also the values of the skewness of dur/dr 12|] which is 


constant and is equal to —3/7. Table 1 reports the comparison between the 
value of the skewness 

^ {{dur/drf) 

’ {{dur/drffP 

calculated with the proposed expression of K and those obtained by the sev¬ 
eral authors with direct numerical simulation of the Navier-Stokes equations 
(DNS), and Large-eddy simulations (LES). It results that the maximum ab¬ 
solute difference between the proposed value and the other ones results to 
be less than 10 %. Other comparisons between the proposed closure for¬ 
mulas and the results of the literature can be found in the previous works 


Q 


,y,ii4|. 
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Table 1: Comparison of the results: Skewness of dur/dr at diverse Taylor-scale Reynolds 
number Rt following different authors. 


Reference 

Simulation 

Rt 

HsiO) 

Present result 

- 

- 

-3/7 = -0.428... 

m 

DNS 

45 

-0.47 

129] 

DNS 

64 

-0.40 

[8J 

DNS 

202 

-0.44 


LES 

00 

-0.40 

LU 

LES 

< 71 

-0.40 

[la 

LES 

720 

-0.42 


Equations (1H|) coincide with those obtained in Refs. [l2, 13| through the 
properties of the motion of the hnite-scale Lyapunov basis and of the frame 
invariance of K and G. This shows that the hypotheses adopted in Refs. 
[ 12 , 13 1 agree with the assumptions of the present analysis. 

The main asset of Eqs. (H4l) with respect to the other models is that Eqs. 


are not based on phenomenologica 


the existence of the eddy viscosity [22 


assumptions, such as for instance, 


23 


25 


15 


3l| . but are obtained 


through theoretical considerations regarding the statistical independence of 
^ from u, and the Liouville theorem. Thanks to their theoretical foundation, 
Eqs. do not incorporate free model parameters or empirical constants 
which have to be identihed. 

Refs. [ 12 , 13 1 show that such formulas describe adequately the energy 
cascade mechanisms. Specihcally, K reproduces the phenomenon of kinetic 
energy cascade in line with the Kolmogorov law, and G describes the thermal 


19 































energy cascade according to the theoretical argumentation presented in the 


Refs. 


26| . to experimental results 
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, 12^ , and to numerical data 


3, 


16|. 


K and G arise from inertia forces and convective term, thus these do not 
depend on v and y. Specifically, K and G are only indirectly related to v and 
X through / and fg whose time evolutions depend on diffusivities by means 
of Eqs. (E]) and ([3]). 

Finally, we remark the limits of the proposed equations. Such limits arise 
from the hypotheses under which Eqs. fl44p are derived: Eqs. 044 p hold only 
in regime of fully developed turbulence where the flow statistical properties 
exhibit homogeneity and isotropy. Otherwise, during the transition through 
the intermediate stages of the turbulence, in decaying turbulence, or in more 
complex situations of developed turbulence with boundary conditions, for 
instance in the presence of wall, Eqs. fl44p can not be applied. 


Remarks. Observe that, thanks to the Navier-Stokes bifurcations, the 
proposed closure formulas modify significantly the mathematical structure of 


Eqs. (j2]) and 
In Refs. 


(H,- 

18 


19| . the authors, studying the non-closed von Karman- 


Howarth equation by group theoretical methods, suggest solutions to the 
closure problem of isotropic turbulence, especially for what concerns the de¬ 
cay of the turbulence. Their two works were not based on the Navier-Stokes 
bifurcations, and their proposed closure formulas exhibit symmetries. Here, 
Eqs. may not present such symmetries, and this is due to the Navier- 
Stokes bifurcations which determine the continuous trajectories divergence 
and thus the possible absence of these symmetries. 
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7. Main properties of the closed equations solutions 

This section analyzes some of the properties of the closed Corrsin and von 
Karman-Howarth eqnations, with particnlar reference to the evolntion times 
of the developed kinetic energy and temperatnre spectra. In particnlar, we 
will show that these spectra reach their developed shape in finite times which 
depend on the the initial condition and on the classical maximal Lyapnnov 
exponent A. 

The exponent A is linked to A throngh 


A = lim A(r) = A(0) 

r-^O 


(45) 


and can be calcnlated with Eq. (I33|) considering that, near the origin / and 
fo, behave like 

2 


fir) = 1 - 


1 / r 


2 V A 


+ ..., 


(46) 


Hence 


feir) = 1 - 


m = 


A. 


+ .. 


ujt) 

Ar(^) 


(47) 


Now, with reference to von Karman-Howarth and Corrsin eqnations, u and 
9 decrease with t according to 


1 1 ilic lll- 

mlai 


du^ 

lOiyu 

dt 

XI 

de^ 

12x9‘ 

dt 



( 48 ) 
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whereas the correlation scales change following these equations 
dXj’ 


dt 


u 

= - 2 +" 




dXe _ 6 


( 49 ) 


where Eqs. fl48|) and fl49|) are, respectively, the equations for the coefficients 
of the powers and of von Karman-Howarth and Corrsin equations. 
Equations fl49)) include the terms with /^^(O) and /g^(0) whose determina¬ 
tion requires the equations for the coefficients of the power of r higher than 
2. Therefore, we qualitatively discuss the variation laws of the correlation 
scales and of the classical Lyapunov exponent. The first terms at the R.H.S. 
of Eqs. fj49l) are responsible for the energy cascade, whereas the other ones 
are due to viscosity and thermal diffusivity. For what concerns the energy 
cascade, it tends to reduce the correlations scales, and if this is sufficiently 
stronger than viscosity and thermal diffusivity effects, then dXx/dt < 0 and 
dXg/dt < 0. On the contrary, the diffusivities tend to increase the correlation 
scales, thus /^^(O) and //^(O) are expected to be positive quantities. 

For sake of our convenience, the condition z/=0, x=0 is first analyzed. In 
this case, u and 6 are both constants, whereas At, Xg and A vary with t. 
From Eqs. fl49|) and fITT)) . At and Xg are proportional with each other, and 
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Figure 2: Taylor and Corrsin scales, and classical Lyapunov exponent in function of the 
dimensionless time. 


vary linearly with the time according to 

AT(t) ^ Xejt) _ 

A'r(O) A0(O) 2’ 

^ (50) 

A(0) l-r/2’ 


r = A(0)t 
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where r is the dimensionless time. Therefore, for z/=0, fhe energy 

cascade described by Eqs. fHT|) determines the decreasing of the Corrsin and 
Taylor scales until to r —)■ 2, where both the spectra are considered to be 
fully developed, and —)■ 0, Ag —)■ 0 and A —)■ cx) (see solid lines of Fig. 

[2]). That is, the correlation functions exhibit developed shapes in a finite 
time whose value depends on the initial condition A(0). The correlations 
scales are decreasing functions of r whereas u and 9 are both constants, and 
this means that mechanical and thermal energies are transfered from large 
to small scales. 

For u >0, X >0, du/dr <0 and dO/dr <0, therefore / and fe are consid¬ 
ered to be fully developed when dXx/dt = 0 and dXg/dt = 0, respectively. 
These cases are qualitatively represented by the dashed lines for different 
values of Rt and Pe, where Rt = Xtu/p, Re = PtRt and Pr = vjx are, 
respectively, the Reynolds number and the Peclet number, both referred to 
the Taylor scale, and the Prandtl number. If the initial values of Xt and 
Xq are relatively high, the dissipation effects are quite small in comparison 
with those of the convective terms. Accordingly, the phenomenon of energy 
cascade is initially stronger than the diffusivities effects, and this determines 
that Xt and Xq diminish exhibiting about the same trend of the case v = x=0. 
Following Fqs. fl48|) and fl49|) . the interval of r G (0, 2) can be splitted in two 
regions for / and fg. The hrst ones occur where dXr/dr <0, dXg/dr <0 until 
to certain values of Ti <2 , T 2 <2, in which dXT/dtiri) = 0, dXg/dt{T 2 ) = 0 
(dashed lines) where, in general, ti ^ T 2 - In this last situation, the mecha¬ 
nism of energy cascade is balanced by viscosity and thermal diffusivity, and 
the turbulent spectra can be considered to be fully developed. In any case. 
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this happens in finite times r <2 for the two spectra. The Lyapunov expo¬ 
nent initially coincides about with that calculated for v =0, then reaches its 
maximum for ri < 2 and thereafter diminishes due to the viscosity effects. 
When A reaches its maximum, dKjdr =0, it is reasonable that chaos and 
mixing achieve their maximum level, and the spectra are there considered to 
be fully developed. 

Thereafter, we observe the region where dA/dr <0. There, due to the 
smaller values of Xt and Ag, the dissipation is stronger than the energy cas¬ 
cade, and Xt and Xe tend to rise according to Eq. fH9|) . This region, which 
occurs immediately after the fully developed condition {dA/dr =0), corre¬ 
sponds to the decaying turbulence regime. 

It is worth to remark that the proposed closure equations flUD are ex¬ 
pected to be valid in the region where dA/dr >0, where the effects of the 
Navier-Stokes bifurcations generate the regime of fully developed turbulence. 
On the contrary, for decaying turbulence, dA/dr <0, after a certain time, say 
r’*' > Ti, the regime of decaying turbulence may not correspond to the fully 
developed turbulence, and Eqs. fHT|) are not defined. 


In the case of complete self-similarity, the several solutions of Eqs. ([2])- 
03 ]) are congruent with each other by proper scale factors depending on only 
one of the variables. Thus, the von Karman-Howarth and Corrsin equations 
are reduced to be ordinary differential equations, and this happens when 

tfl, 


dXll 


dt V 


— = const. 


—^ = Rt = const, = const 
V Xt 


(51) 


On the other hand, the proposed closure formulas can not be brought to 
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Eqs. dSU), thus Eqs. (m|) do not give a complete self-preservation. Never¬ 
theless, in the cases in which the dimensionless quantities of Eq. fl5T]) exhibit 
relatively slow variations with respect to the correlations, the solutions can 
be considered to be self-preserved only in Erst approximation. In this case, 
the correlations read as 


/ = / 





9{t) u{t) 

Next, from Eq. fl48ll we have 

\g{t) ^ 

Xrit) \ 5 Pr 

and is proportional to ^ according to Eq. fH9|) 

0 ) 7 , 

/'>'((, 0) 3 


( 52 ) 


( 53 ) 


( 54 ) 


Hence, the correlations are considered to be fully developed when t* = r-i = 
T 2 . There, f^^{0,t*) is linked to Xrit*) through Eq. (14^ with dXT/dt=0 14\ 




7 X^{t* 


Rt{P 


+ 5 > 0 


(55) 


We conclude this section by observing that the times of developing of / 
and fg are hnite quantities which depend on the initial conditions. 
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8. Conclusions 


In order to obtain the closure formulas for von Karman-Howarth and 
Corrsin equations, this analysis represents the fluid motion in the Lagrangian 
form. As the separation vector varies much more rapidly than the velocity 
field, ^ and (u, '&) are assumed to be statistically independent. This con¬ 
jecture and the adoption of the Liouville theorem lead to the closure of von 
Karman-Howarth and Corrsin equations. The equations here obtained coin¬ 
cide with those of the previous study [l2|, ll3|, showing that the approach of 
Refs. 12|, [iSj], dealing with the frame invariance of K and G and the prop¬ 
erties of the finite-scale Lyapunov basis, is in agreement with the present 
analysis. This corroborates the previous works, providing a demonstration 
of the closure formulas completely different with respect to the previous ar¬ 
ticles. 

Thereafter, some of the properties of these equations are studied. In 
particular, we show that the condition of developed spectra (correlations) 
in isotropic homogeneous turbulence is reached in hnite times whose values 
depend on the initial conditions. 
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